Linear Growth of Structure in the Symmetron Model 
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In the symmetron mechanism, the fifth force mediated by a coupled scalar field (the symmetron) 
is suppressed in high-density regions due to the restoration of symmetry in the symmetron poten- 
tial. In this paper we study the background cosmology and large scale structure formation in the 
linear perturbation regime of the symmetron model. Analytic solutions to the symmetron in the 
cosmological background are found, which agree well with numerical results. We discuss the effect 
of the symmetron perturbation on the growth of matter perturbation, in particular the implications 
of the brief period of tachyonic instability caused by the negative mass-squared of the symmetron 
during symmetry breaking. Our analysis and numerical results show that this instability has only 
very small effects on the growth of structures on sub-horizon scales, and even at horizon scales its 
influence is not as drastic as naively expected. The symmetron fifth force in the non-tachyonic 
regime does affect the formation of structure in a nontrivial way which could be cosmologically 
observable. 



I. INTRODUCTION 

Scalar fields coupled to matter are generic predictions 
of many theories of high energy physics. In recent years, 
this idea has attracted a lot of attention in the context of 
dark energy [1], which is believed to be a scalar field or 
one of its variants [2-6]. However, it is well known that 
if a scalar field couples to matter or curvature then a 
scalar fifth force and a modification to the gravitational 
law could result. Such new physics has been strongly 
constrained by local gravity experiments and solar sys- 
tem tests, so that fifth forces and modifications of gravity 
must been either very short-ranged or very weak, or both. 
If this is true, then the effect of the scalar field is mainly 
to drive the accelerating expansion of the Universe, such 
as in the quintessence model [2]. 

More cosmologically interesting models could be built 
if the fifth force or modification to the standard Ein- 
stein gravity is only weak and/or short ranged where lo- 
cal experiments are performed, such as in our Solar sys- 
tem where matter density is high and gravity is strong, 
but could become strong (of gravitational strength) and 
long-ranged elsewhere. Over the past few years, several 
models have been proposed to realise this, including the 
chameleon model [7-13], the environmentally-dependent 
dilaton model [14, 15], the DGP model [16], the galileon 
model [17-20] and the symmetron model [21, 24, 26]. 

Although all these models predict strong environmen- 
tal dependence of the fifth force, they work in very dif- 
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ferent ways. In the chameleon model, for example, the 
fifth force is suppressed exponentially in high density re- 
gions where the scalar field acquires a heavy mass via its 
coupling to matter. The environmentally dependent dila- 
ton model, on the other hand, drives the scalar field to 
some critical value <p c in high density regions, which cor- 
responds to a vanishing coupling and therefore vanishing 
fifth force (though the mass of the scalar field depends 
on the local matter density as well in this model). 

The symmetron model, which is the topic of this work, 
relies on a similar mechanism to suppress the fifth force 
in high density regions. Here, when the matter density 
is high enough, the effective potential of the scalar field 
has a global minimum at the origin, which corresponds 
to a vanishing coupling to matter. When matter den- 
sity drops below some critical value, the symmetry in 
the effective potential is broken and two local minima de- 
velop and move away from the origin, corresponding to 
a nonzero coupling to matter and thus a non-vanishing 
fifth force. If later the matter density inside a region be- 
comes high again due to the structure formation process 
(such as in galaxies and galaxy clusters), the symmetry 
of the effective potential could be restored and the fifth 
force vanishes again for that region. 

Local experiments could constrain the parameters of 
the symmetron model. Interestingly, in this constrained 
parameter space, there are still models which could devi- 
ate from the standard ACDM cosmology on large scales. 
The cosmological obscrvablcs, in particular those rele- 
vant for the large scale structure formation, can thus pro- 
vide valuable information about the symmetron, such as 
whether it exists, what its observational signatures are, 
how to differentiate it from other models, etc. In this 
work, we shall concentrate on the structure formation of 
the symmetron model in the linear perturbation regime, 
and show how the symmetron can affect the large scale 
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structure of the Universe. 

This paper is arranged as follows: we first briefly 
overview the symmetron model and how the local tests 
constrain its parameters in Sects. II A and II B, followed 
by analytical and numerical study of the background cos- 
mology in the model (Sect. II C) and a short description 
of the tachyonic period, during which the scalar field has 
a negative mass-squared (Sect. II D). Then in Sect. Ill we 
study the general behaviour of the linear perturbations in 
this model, where we show that the fifth force essentially 
enhances gravity within its range (given by the comov- 
ing Compton length of the symmetron field) and keeps 
standard gravity unmodified beyond that range. We also 
show that the tachyonic instability docs not affect the 
structure formation significantly on sub-horizon scales. 
In Sect. IV we give some numerical results of the large 
scale structure formation and show that significant de- 
viations from the ACDM paradigm can be found, which 
make the model cosmological interesting. We summarise 
and conclude in Sect. V. Throughout this paper we shall 
adopt c = h = 1 and the metric convention (—,+,+,+). 



II. SYMMETRON DYNAMICS 
A. The symmetron 

Scalar-tensor theories are characterised by their cou- 
pling to matter and their interaction potential. In such a 
context, the symmetron model was proposed in [21, 24] 
and is described by the action: 
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where <j> is the scalar field and V(4>) its potential, S m = 
S m [vP 1 , <7/ii>] the matter action with the matter fields 
which are minimally coupled to the Jordan frame metric 
Qiiv = A 2 ((j))g^ v ] g^ v is the Einstein frame metric, which 
is used to compute the Ricci scalar R; n = 8irG = m~ 2 
where G is Newton's constant and m v \ the reduced 
Planck mass. 

For the symmetron model, the interaction potential 
and the coupling function are simply chosen such that 



v{4>) = v - \y<t> 2 
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in which Vo is a cosmological constant and has mass di- 
mension 4, p and M have mass dimension one and A is a 
dimensionless model parameter. We shall see that Vo is 
needed to explain the recently observed accelerating ex- 
pansion of the Universe, but being a constant it bears no 
influence on the dynamics of the symmetron field. Note 
that the potential has two nontrivial minima while the 
coupling function is monotonously increasing. 



The field equations are obtained by varying the action 
S with respect to the symmetron field 4>, and we have 



D^=V^)-A^)A 3 (<P)T, 
f = g^f^, 
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where we have defined the Jordan frame energy momen- 
tum tensor that is related to the Einstein frame one 
by T*„ = A 3 {ct))f^ p g pv . Note that we raise and lower the 
indices of using the Einstein frame metric, g^ v . With 
this definition, the field equation for <fi and the Einstein 
equations become 
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in which the symmetron field is governed by an effective 
potential 



V- off (</>;T) = V(0)-A(0)T, 
and the total energy momentum is given by 
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which satisfies the following conservation equation 



(10) 



Note that this implies that for pressureless matter with 
T pl/ = pu p u v with u p a 4- velocity (u^u^ — —1), we have 
—T = p m and p m is conserved independently of <f>: 
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As is the usual practice, we define the effective mass of 
the symmetron field by 
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where we have used the fact that, from the above equa- 
tions, the effective potential can be rewritten as 



(13) 



Hence, in the symmetron model, as long as p m is high 
enough, namely p m > p+ where 



p 2 M 2 
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the minimum of the effective potential is at the origin 
(0 = 0). In contrast, in vacuum, the symmetry is broken 
and the potential has two nonzero minima: 
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As long as </> <§C M the effective coupling to matter reads 

(16) 



a O) = m pl [lnA] !0 w 



leading to the absence of modification of gravity in dense 
environments where the field vanishes. The modification 
of the growth of structure depends on 



a = 2al, 



(17) 



and the order of magnitude of a in vacuum, 

a* = a(0*), (18) 

is crucial for the growth of the large scale structure. For 
values of the density lower than p*, gravity becomes mod- 
ified in a way which could be tested cosmologically, and 
a* characterises the relative strength of the modification. 



It is convenient to simplify this equation by separating 
two regions with different behaviours [21] 



and 



Kff(r>i2) = M 2 (0-^) ! 
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where one assumes that the density vanishes at infinity 
and outside the body. The solutions read 



and 



(r<R) = 4> ln {r) = A^sinh(rJ^-^ 



(r>R) = <f> out (r) = B-e~^ r + ^ (25) 



B. Gravity tests 

In this section, we review some of the gravitational 
properties of the symmetron[21] The symmetron model 
is designed to induce modifications of gravity which could 
be tested in the near future, both gravitationally and cos- 
mologically. Requiring that the energy density at which 
the curvature at the origin of the potential changes sign 
(roughly when gravitation undergoes a transition from 
standard to modified) is close to the current critical en- 
ergy density (which implies that gravity is modified cos- 
mologically recently), we have the estimate 

Moreover, the modification to gravity is detectable only 
if it is comparable to (or bigger than) standard gravity, 
or equivalently the effective coupling ~ C(l), which 
implies 
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M 

where we have defined A 2 
the vacuum mass 

m 2 (^o) = 2p 2 - O 
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H 2 = O{A 2 )H 2 (19) 



and correspondingly set the interaction range of the sym- 
metron to be ~ O (m _1 (0o)). This also determines the 
self coupling [21] 
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It is then crucial to have an estimate for M. This follows 
from the study of solar system tests. 

Let us consider a spherical object of density p and ra- 
dius R. The static solutions for the field profile are ob- 
tained solving 
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Defining the modified Newton potential at the body's 
surface 



P R 2 
M 2 



= 6A 2 $ 



N 



(26) 



and the ratio of the size of the sphere to the range of the 
symmetron interaction 



we find that 
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Notice that a ^> f3 2 as long as p p+. There are two 
types of solutions depending on the values of a. When 
a <C 1 we have 
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The scalar force acting on a test mass outside the sphere 
is 
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(31) 
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when the Newtonian force is -^ £L - = „ p f 2 implying that 
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This implies that the scalar force is not screened. On the 
other hand when 5»1 we have 
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The scalar force is now 
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hence 
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leading to a large screening of the scalar force. 

Phenomenologically, one must impose that in our 
galaxy ^> 1. Imposing at least ac > 10 and upon 
using $ G ~ 10~ 6 , one gets M < 10" 3 m pi or equivalently 



M > io 6 



(37) 



This implies that jf < 10 



The range of the symmctron in vacuum is given by 

1 Mpc (38) 



At" 1 < lO 3 ^" 1 



which corresponds to relevant scales for astrophysics. We 
will come back to this point later. If M « 10~ 3 mp; 
then the scalar field is just about screened by the sun as 
$0 ~ 10~ 6 => a© ps 10. On the other hand, the earth 
is not screened as $® ~ 10~ 9 and a ffi ps 10~ 2 . What 
matters then for solar system tests is the value of the 
field in the galaxy: 



M 



M R G 
mpi tVogR, 



■ exp 



Rg — Rs 
Rg 



lOLG 



(39) 



The most stringent constraint in the solar system is the 
Cassini bound I7— 1| < 10 -5 on the Eddington parameter 
7 [22]. In the Einstein frame, the metric is expressed as 
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while the one that particles follow is the Jordan frame 
metric 



5oo 
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9ij = (l-2 7 $j)%, 
from which we have [21] 
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Denoting by the moment when the symmetron poten- 
tial becomes unstable at the origin, we have 



p?M 2 
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(45) 



Fixing for instance ac ~ 600 and M ss 10 _4 mp; while 
using for the galaxy /3 ~ 4 x 10~ 2 (1 + z*)^ we get a 
bound on admissible z+ < 20 to evade the Cassini bound 
|7- 1| < 10~ 5 . 



C. Cosmological evolution 

We are interested in the symmetron evolution during 
the matter dominated era, when the symmetron potential 
is negligible compared to the matter density p al . This 
can be seen by evaluating the value of the potential at 
its absolute minimum in vacuum as the potential at 
the effective minimum is monotonously decreasing with 

Pm 

V(M - -£ 

which is always < 10 -6 HQm 2 pi . Moreover, as 4>+/M < 
10 -3 , the coupling function A(<j>) is also approximately 
equal to one. The Friedmann equation is then to a high 
precision 



H 2 = 



3m 2 Pl 



(46) 



in the matter dominated era. 

A full picture of the time evolution of the reduced sym- 
metron field ip = 4>/4>-k is shown in Fig. 1, and we can see 
that after the symmetron breaking at x{= t/t^) = 1 the 
symmetron deviates from zero, and finally starts to os- 
cillate around the moving minimum of the effective po- 
tential (V c ff). The detailed evolution of <p is of course 
more complicated, but still understandable analytically 
in certain limits, as we shall show now. 

Let us denote by the instant when p m = p+ and a* 
the corresponding scale factor. The scalar field equation 
of motion reads now 
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which can be cast in a dimensionless fashion by defining 
k = jit* in addition to the tp and x defined above, where 
we have the limit 
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The dimensionless version of the symmetron equation of 
motion is then given by 

ip"{x) + -i P '{x) + n 2 (\-l]Lp{x)+K 2 ip{xf =0, (49) 
X \x z / 

where ' = d/dx. 

Before the curvature at the origin changes its sign, the 
symmetron oscillates about the origin <f> = 0, which is 
then the minimum of the effective potential. Assuming 
that the initial amplitude is small, the symmetron equa- 
tion of motion simplifies to 



x 2 ip" + 2x(p' — k 2 (x 2 — l)ip ps 



The solutions read 



(p(x) 



[aiI iL {nx) + biK u (nx)} 



(50) 



(51) 



5 




1,10 



FIG. 1. The time evolution of the background symmetron field tp = </>/</>*, in which the parameters are chosen as k = 4 x 10 , 
Zi, = z(t*) = 3 and a* = 1. Note that immediately after the symmetry breaking, i.e., x = 1, tp grows exponentially but then 
quickly begins to oscillate. The oscillation is around the minimum of the effective symmetron potential V e s{4>). Note also that 
—tp is also a possible branch of solution but it gives nothing new. 



where t = V4 ^~ 1 



a\ and 62 are constants of integra- 



tion which could be determined by the initial conditions. 
As long as 1 < 1 we can expand 



m = mi j 
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The upper left panel of Fig. 2 shows a numerical example 
of the evolution of tp with respect to x for x < 1 (t < t*), 
from which we can see the oscillation around the global 
minimum of V^g (</>): (j> = 0. 

After the change of curvature, the field rolls away from 
the origin and lags behind the minimum of the effective 
potential. The minima of the effective potential are at 
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Before reaching one of the minima and oscillating around 
it, the symmetron field will first linger around the origin 
before following the inflection point where the curvature 
of the effective potential vanishes and then settling down 
to that new minimum. At the inflection point the field 
is close to 
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where 



0, and in deriving this equation we have 



used the facts that p m a 3 = p+a 3 (conservation of mat- 
ter) and that in the matter dominated era a oc t 2 ^ 3 . At 



this time, the first derivative of the symmetron effective 
potential becomes 
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The symmetron equation of motion becomes then 



tp + -tp 
x 



2k 2 
3\/3 



3/2 



(55) 



the solutions of which arc 
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where 02 and hi are again constants of integration which 
depend on the initial conditions. The upper right panel 
of Fig. 2 shows the comparison of the numerical solu- 
tion with the above analytic approximation (51) , where 
we can see that for small x just over 1 the agreement 
is very good. In the lower left panel of the same figure, 
one shows the short period of time when the symmetron 
tracks the inflection point (50) before reaching the mov- 
ing minimum. 

In a third phase, the symmetron field catches up with 
the minimum of its effective potential. Let us define tp = 
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FIG. 2. The numerical results and analytic approximations for the background evolution of cp = 4>/4>* m different stages, for 
which the model parameters are chosen as k = 4 x 10 3 , z* = z(t*) = 3 and a* = 1. The initial condition for the symmetron 
is chosen as tp(z eq ) — 1CF 2 . Upper left panel: the numerical oscillation solution to ip before the symmetry breaking, i.e., when 
x < 1 or equally t < t*. Upper right panel: the numerical (red curve) and analytic (green curve) solutions to ip immediately 
after the symmetry breaking, when the scalar field lags behind the minimum of its effective potential V e g. Lower left panel: 
The analytical (green curve) and numerical (red curve) solutions when the symmetron tracks the inflection point. Lower right 
panel: the numerical (red curve) and analytic (green curve) solutions to <p after the scalar field settles down to the minimum 
of V e B and starts to oscillates about it. 
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Around the minimum of the effective potential, ip + (t) 
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one can expand to first order ip = ip 
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This is a forced oscillation where the forcing term is neg- 
ligible for x 3> 1 and the mass-squared term satisfies 

m 2 = 2fi 2 ^1 — j^J The time variation of the mass com- 
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pared to the period of the oscillation is given by: 
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which is small when <C 1. a condition satisfied when 
|t — 1 3> k~ 3 , i.e., large enough x. The solution for large 
enough a: is then given by 



t 2 -a 



a 3 cos O(0 + b 3 sin Cl(t)] , (59) 



where 



Q(t) = V2k iJ^-l + taxi' 



and 03 and 63 are constants of integration which depend 
on the initial conditions. 

A better solution can be obtained by taking into ac- 
count the forcing term whose characteristic time is r = 



£± 

<i>+ 



The forcing term is slowly varying when 



— < 1, which reads 
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and is indeed small when the period is slowly varying. In 
this better approximation is given by 



,(1) 



(60) 



where the forcing term rapidly becomes negligible. The 
lower right panel of Fig. 2 shows the analytic and numer- 
ical solutions for this stage; note the drift in the period 
of the oscillation due to the fact that close to x = 1 the 
ratio fn/m 2 is not negligible. 



D. Tachyonic instability 

The mass-squared of the symmetron field is briefly neg- 
ative when the symmetron lags behind the minimum of 
its effective potential. Defining 

in which tf is the time when the symmetron field set- 
tles down to the effective potential minimum, one could 
estimate At/t+. In the case for k = 10 3 , a good approxi- 
mation for the time spent in the tachyonic regime can be 
numerically fitted as 

w -2 x 1(T 4 log 2 (ip(t eq )) - 9.7 x 1(T 3 \og(<p(t eq )) 
+2 x 10 -2 . 



As it stands, this tachyonic period is extremely short (see 
the discussion below) . We will see shortly that the effect 
of the negative mass-squared is only relevant on large 
enough scales for cosmological perturbations, and even 
for those large scales it is quite insignificant. 

III. LINEAR PERTURBATIONS 

A. Growth of structure 

Wc will now be interested in the growth of linear per- 
turbations in the matter dominated era, and for that we 
work in the Newtonian gauge where the perturbed metric 
in the absence of any anisotropic stress is given by 

ds 2 = -(1 + 2$ N )dt 2 + a 2 (t)(l - 2§ N )dx i dx i (61) 

where $at represents the Newtonian potential. For sim- 
plicity let us assume that matter comprises a single fluid 
of pressureless particles with the energy momentum ten- 
sor v = Apu^u v , p being the conserved energy density 
and w M the four- velocity of the matter particles. In gen- 
eral the conservation equation reads 



^ + 3hp = 
ar 



(62) 



where t is the proper time along the particle trajectories. 
The local Hubble expansion rate is defined as 3h = V^u^. 
The perturbed conservation equation is then 
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The Euler equation in terms of the divergence = V • u 
of the velocity field becomes here 
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which is modified by the presence of the symmetron field 
as indicated by the last term on the left-hand side. 

Similarly, the modified Poisson equation now involves 
the perturbation of the Einstein frame matter (and scalar 
field) density p = A((f>)p and reads 
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which becomes in Fourier space when A{<f) « 1 
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The scalar field equation of motion is expressed as 

= -2® N Vi ff {4>) + 4$W - £td (67) 

From the above equations one can derive a second order 
differential equation for 8 = with the coupling to the 
scalar field taken into account: 
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where have defined 8 = — — . 

r m p i 



B. Adiabatic approximation 



the symmetron lags behind the minimum, we find that 
gravity is modified according to the comoving Compton 
radius am with an amplitude depending on a. Structures 
on scales outside the Compton radius grow as in GR 



We have seen above that the symmetron field lags be- 
hind the minimum of its effective potential for a very 
brief period when the mass-squared m 2 (</>) < 0. In such 
a tachyonic phase, the perturbation 5(f) grows exponen- 
tially for modes such that |y + m 2 {4>) < [26] as can be 
seen from Eq. (67). At the time of symmetry breaking, 
the effective mass-squared vanishes and then decreases 
to a fraction of — before increasing to its value at the 
effective potential minimum. Clearly larger-scale modes, 
for which k is smaller, are easily in the tachyonic regime; 
as \x is much larger than the Hubble expansion rate today, 
some of the tachyonic modes could indeed be sub-horizon 
(big k) and this could possibly influence the growth of 
large-scale structure on sub-horizon scales. 

Outside the tachyonic regime, the evolution equation 
of the symmetron perturbation is sourced by the matter 
perturbation. Assuming that the symmetron field tracks 
the effective potential minimum, i.e., for non-tachyonic 
modes and neglecting the short period during which the 
symmetron field lags behind the minimum, its equation 
of motion in the sub-horizon limit becomes 



(69) 



where the time derivatives are much smaller than the 
spatial derivatives in the sub-horizon limit and thus are 
neglected. From this equation the symmetron field per- 
turbation can be solved as 
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which shows that the symmetron field perturbations 
tracks that of the matter density. This in turn implies 
that matter perturbation well within the horizon grows 
according to 

. . S k 2 

6 + (2H + p<p)5 - -H 2 S — =-/3 60 
2 a z 



or equivalcntly 



5 + 2H6 - -H 2 



l + a{4>)- 



1 



a 2 m 2 ((f)) 

— p — 



8 = 0. (71) 



Therefore, in this adiabatic approximation valid for non- 
tachyonic modes and neglecting the interval of time when 



while those on scales inside the Compton radius have a 
modified growth due to the renormalised Newton con- 
stant 



G c ff = G N (1 + a). 



(72) 



We will see in the following that this result is hardly 
modified by the tachyonic instablity. 



C. Tachyonic instability 



As mentioned earlier, if the mass-squared of the sym- 
metron field m 2 (0) becomes negative (see Fig. 3), then 
the perturbation will undergo an unstable growth which 
could be problematic in some cases. We have also seen 
that such tachyonic instability problem is most likely to 
plague the large-scale (small k) modes. The purpose of 
this subsection is to assess how big the impact it could 
have on the growth of matter perturbations. 



0,10- 



0,04- 




1,015 



FIG. 3. The evolution of the mass-squared of the symmetron, 
from which we can see that for a short period immediately 
after the symmetry breaking (x = 1) the mass-squared be- 
comes negative and the symmetron field and its perturbation 
experience tachyonic instabilities. 
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We are interested in modes which enter the horizon af- 
ter matter-radiation equality and before the tachyonic in- 
stability happens. Normalising a* = 1, this corresponds 
to 

! >1 (73) 



and 



k 



l + z e 



< 60 



(74) 



1 + z+ 

The tachyonic modes can be conveniently studied using 
the reduced symmetron perturbation 

x ^ 
dip = — 



At the horizon crossing the symmetron perturbations are 
taken to vanish, as during inflation the mass of the sym- 
metron is much larger than the Hubble expansion rate. 
This implies that the initial conditions for the symmetron 
perturbation are 

Scf>t k « 0, 8<j> tk « (77) 
where tk is the horizon-entry time for the mode k 



tk 



TJkHl 



(78) 



Using the fact that the tachyonic growth starts at x = 1, 
we have 



and the parameters 



< 10-3 
M ~ 



and 



2k 



3#* 

Keeping terms in ui 2 (or f3), neglecting terms in lo a (or 
/3 2 ), and considering that the variations in both ip and 
Sip are very rapid in the tachyonic period, the growth 
equation simplifies to 

4 2 2ui 2 

S" + —5' - ^S = -—^foStp" + 2ip'Sip' + tp"5<p) 



3x 



3a; 2 



in the sub-horizon limit for the tachyonic modes. Notice 
that in the absence of the rapid variation of Sip due to the 
tachyonic instability, this equation reduces to the growth 
equation in GR. 

Defining S = x»5, the growth equation becomes 

S" -^8* -%(<pSip" + 2p'8p' + p"5 v ), (75) 

VX !; k 

the solutions of which are 



Oi\X 3 + 



x s 



6^ 

















-4) 


Jx k 


\X3 


S3 / 



{<p6<p)"(s)ds, 



where u\ and f3\ are integration constants, and Xk is the 
value of x when the fc-mode under consideration enters 
the horizon. As a result, 



5(x) = a,\x~i 



Integrating by parts 




(ipSip)" (s)ds. 



r (pSip)"(s)ds 

X S3 / 



£1 _ £1 {ip5ip y {s) 

X S3 




(ipS<p)'(s)ds. (76) 



4 S3 

3 x 

4 S3 
3 x 



1 X3 

3 si 

1X3 

3 si 



(ipSip)' (s)ds 
(pSip)' (s)ds 



(79) 



The variation of p>Sip is much faster than the other terms 
in this integral, so that the terms in the brackets could 
be absorbed into the derivative with respect to x, and 
the integrand becomes a total derivative. This implies 
that the growth factor behaves like 



S(x) 



if x < 1 



(80) 



As structures grow anomalously due to the tachyonic in- 
stability, we define the deviation from the normal growth 
as 



A(a:) = S{x) ~ S h 



(81) 



When x > 1, this becomes (remember that we have nor- 
malised a* = 1) 



2u? 



■ipSip(x) 



9 a 2 H 2 
'2~k?~' 



(82) 



Our numerical and analytic solutions to A are shown in 
the lower left panel of Fig. 4, and the agreement is very 
good. To see how big an impact the tachyonic instability 
could have on the growth of matter density perturbation, 
let us consider the quantity 



A 

Sgr 



9 a 2 H 2 



2k 2 S 



GR 



Scj) 
'rripi 



"tU" ) ^ (83) 



where in the second step we have used the Poisson equa- 
tion Eq. (66) in GR, i.e., with the Sip term removed. To 
have a rough idea about the magnitude of this quantity, 
we need the following ingredients: 
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0,03- 



0,02- 



-0,02- 




2, x 10 11 



l.x 10 11 



-l,x 10"" 



FIG. 4. The numerical results and analytic approximations for the perturbation evolution in the symmetron model, for which 
the model parameters are chosen as k = 4 x 10 3 , jz* = a(t*) = 3 and a* = 1. The initial condition is chosen as ip(z eq ) = 1CF 2 . 
Upper left panel: analytic (green curve) and numerical (red curve) solutions to the symmetron perturbation Sep = o~4>/4>*\ note 
the excellent agreement between the two. Upper right panel: the time evolution of (p (green curve) compared to that of 8(p; 
note that \6<p\ is always smaller than ip, and does not blow up exponentially. Lower left panel: the analytic (green curve) and 
numerical (red curve) solutions to the A defined in the text; note the excellent agreement between the two. Lower right panel: 
the ratio A/S which is much less than unity, showing that the correction to the growth of matter perturbation due to the 
tachyonic instability is negligible. 



1 . During the brief period of tachyonic instability, not 
only 5<p but also the symmetron field <f> itself grows 
exponentially. In particular, the mass-squared for 
4>, m 2 (</>), is more negative than for <$</> because the 
latter is k 2 /a 2 +m 2 (4>) > m 2 ((j)) for all fc-modes. As 
a result, although the growth of 5<f> is unstable, it 
is still slower than that of <f> and we expect 8<f>/ <j) < 
0(1), as can been seen from the upper right panel 
of Fig. 4. 

2. We have seen above that for the symmetron model 
to evade the local experimental constraints we must 



have (f>/M < <p*/M ~ 0(1O~ 3 ), as well as M/m p \ < 
0(1O~ 3 ). As a result, <£/m pl < 0(1O~ 6 ). 

3. k/aH 1 for the sub- horizon modes and k/aH ~ 
0(1) for near- horizon modes. 

4. $a?,gr is the gravitational potential in general rel- 
ativity, which should be much less than unity, and 
we could use ~ 0(1O~ 7 — 10~ 3 ) for conservative 
estimates. 

5. o*~0(l). 



11 



For sub- horizon modes, taking k/aH ~ 50 and $jv,gr ~ 
10 -4 , we find that |A/<5 GR | ~ O(10" 8 ) (the exact value, 
of course, depends on the values of (f>, M and 5(f)), which 
is consistent with the numerical result shown in the 
lower right panel of Fig. 4. For near horizon modes, 
on the other hand, k/aH ~ 0(1) and thus we have 
|A/5qr| ~ 0(1O~ 6 /^JV,Gr), i-c, the very large scale 
growth of matter perturbation could be modified sig- 
nificantly and leave potentially observable signatures 1 if 
I^n.grI ~ 10~ 6 . In the present paper, however, we will 
be primarily interested in the sub-horizon scales and ne- 
glecting the tachyonic instability is a reasonably good 
approximation. 

In fact, a very accurate approximation for Sip can be 
obtained by noticing that <p satisfies 



3Hcf) 



~H 2 + m 2 (0) 



while 5d> is a solution of 



5<j)+2,H8cf>+[ ^ + m 2 (0) ) 5<t> 



-2§ N V\4>)+A$ N $-^8 
(85) 

As H 2 is small compared to m 2 ((f)), 8(f) and (f> satisfy the 
same equation provided ^ -C m 2 ((f)) and one can ne- 
glect the right hand side of the Klein-Gordon equation. 
As m 2 ((f)) = 0(n 2 ~ 10 6 iJg) and the tachyonic modes 
are smaller than k < 60ii*, the first condition is satis- 
fied. Moreover, the right-hand side is negligible provided 
8 <C Sip. In this case 8(f> and (f> satisfy the same equa- 
tion provided 5 is small compared to As a result we 

have [23] 



Sip W -r-Slfi, 



(86) 



very brief, thanks to the exponential growth of which 
quickly settles down into the non-tachyonic regime; as a 
result the blow-up of cither (f> or S(f> simply cannot ac- 
cumulate enough momentum to have big impacts on the 
structure formation. 

The study of the structure formation of the symmetron 
model in nonlinear regime is beyond the scope of present 
paper and numerical simulations have been performed in 
[26] . Be aware that an accurate study of the tachyonic 
behaviour of the model will involve a very big simulation 
box (to include the near horizon scales) and the full tem- 
poral evolution of 5(f), which is not achievable by known 
A-body simulations at present. 



IV. STRUCTURE FORMATION IN THE 
ADIABATIC APPROXIMATION 

Most of the sub-horizon scales are in the non-tachyonic 
regime while the symmetron field oscillates around the 
(moving) minimum of the effective potential. This means 
that the adiabatic approximation, where the scalar field 
(in background cosmology) always follows the minimum, 
provides a good description of perturbation theory. This 
has been verified by numerical simulations [26]. In this 
section we shall study the linear structure formation 
quantities in the symmetron model under the adiabatic 
approximation. 



A. The Adiabatic approximation 

In the adiabatic approximation, we neglect the effect of 
the short time period when the symmetron perturbation 
becomes tachyonic which, as we have seen, is justified. 
In this case wc have 



This approximation is numerically very accurate for all 
H* ^ k i5 60 if* as can be seen in Fig. 4. 

Why does the exponential blow-up oi 8(f) have small 
effects on the growth of matter perturbation? First, 
the rapid growth of 5(f) is associated with an even faster 
growth of (f> and so the symmetron perturbation does not 
really enter the nonlinear regime (see upper right panel of 
Fig. 4). In the background cosmology, the blow-up of (f> 
has not spoiled the matter dominated era, and certainly 
we should not expect an even slower perturbation growth 
to be too much problematic in the linear regime. Second, 
as mentioned many times above, the tachyonic period is 



1 Note that even on horizon scales where the tachyonic instability 
is supposed to be more severe, we will not have an unbounded 
blow-up of the matter perturbation as naively suggested by the 
equation of motion for the symmetron perturbation. Note also 
that near the horizon the approximations used to derive the re- 
sults in this subsection might be inaccurate, but we expect them 
to provide a reasonable order of magnitude estimate. 



0min(p) = <f>*\jl - — 9 (P* - P) 



(87) 



o B i,(p) = a*(l-^|%-p), (88) 



where 9(x) is the Hcavisidc function, and 



m mi„(p) 




(89) 



The symmetron model is characterised by the dimension- 
less parameters z*, A 2 and a* determining the redshift 
at the transition 

P* 

PmO 



(1 + ^) 3 



(90) 



In terms of these parameters, [J, 2 = 3A 2 ri m o-ffo (1 + z *) 3 
where f2 m offo = K l/ 5 mo/3. We already know that local 
tests of gravity require A 2 > 10 6 . For z+ > this implies 
that /I 2 > 10 6 7Jq 3> H§ and therefore the symmetron 
mass is always much larger than the Hubble rate. 
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We define 4>(t) to be the value of the symmetron held 
in the cosmological background. The symmetron equa- 
tion of motion and Einstein equation in the cosmological 
background reduce to 



K 2 

H 2 = -y [Pm{a) + /Orad(a) + V + p<j\ 



(91) 
(92) 



P4> = 4^ pm{ - a) - \/~^ + + 5^-( 93 ) 

where p m (a) oc a~ 3 and p xa d(a) oc a~ A . Because p 2 ^> 
H 2 today, one finds that in the cosmological background 
P~4> <C Pm and in the adiabatic approximation the cosmo- 
logical symmetron <p tracks the minimum of the effective 
potential, which is a cosmological attractor: 



4> ~ 0min(/O m ) = 4>*\ 1 



Pn 



<(p*-p m ). (94) 



In terms of the redshift z, this becomes 



'1- 



l + z 

1 + 2* 



(z* - z), 



(95) 



Since p^ is negligible, the background cosmological be- 
haviour of the symmetron model is just the unmodified 
ACDM model with a cosmological constant A = k 2 Vo: 

H 2 ~ ~ [Prad(a) + Pm(a)} + y 

We define m((f>) = fh(z) and a(<f>) = a(z) in the cosmo- 
logical background and then we have: 



m 2 (z) 



a(z) 



3A 2 n m0 H 2 [(1 + zf - (1 + z+Y] z > z*, 

6A 2 n m0 H 2 [(1 + z*Y - (1 + zY] z < z*. 




( 1+* 



z > z*, 



z < z+ 



(97) 



This allows one to study the growth of structure in the 
adiabatic approximation. 



B. Linear Structure Formation 

We assume that high density localised clumps of mat- 
ter, e.g. galaxies, have p 3> p* and therefore are ef- 
fectively screened from the symmetron force. We define 
5 g (x, t) as the relative large scale coarse grained density 
perturbation in the galaxies, and <5 m (x, t) to be the large 
scale (linear) density perturbation in the smooth cold 
dark matter. We work to linear order in <5 g <C 1 and 
5 m <C 1. We define £ g (k, t) and <5 m (k, t) to be the Fourier 
transforms of the relative density perturbations. 



Using the adiabatic approximation for the background 
evolution of the symmetron field, and taking p = — ln(l + 
z), we have (see e.g. [25]): 



+ 



a{z) 



where p = d/dp, and 

relpm(a) 



H 2 



1 + (1 + Z)- I fc- 2 7T1 2 (Z) 



Ml + z?H 2 
H 2 



Galaxies are often used as observational tracers of the 
linear dark matter perturbation. The latter feels both 
gravity and the fifth force due to the symmetron coupling 
whilst the former only evolves under gravity. This is very 
similar to the chameleon and environmentally dependent 
dilaton scenarios [7, 8, 10-13]. Working in the Fourier 
space, we have solved the relevant perturbation equations 
numerically to find the relevant quantities (for simplicity, 
in what follows we shall neglect the tilde of 6 g and S m ) 



• The growth rate / ga i = d(ln<5 g )/dlna. 

s , measured by weak gravita- 



• The slip function, E K 
tional lensing. 



• The slip function measured by the integrated 
Sachs- Wolfe (ISW) effect. 

• The indicator of modified gravity Eq. 

Let us first recall their definitions. In a modified gravity 
setting, there are two metric potentials $ and "f: 

ds 2 = a 2 { V ) [-(1 + 2*)d?? 2 + (1 - 2<5>)dx 2 ] 

In the symmetron model here, the absence of anisotropic 
stress implies 



N 



as in Eq. (Gl). Weak lensing is sensitive to $+\E f while the 
ISW effect is proportional to $ + vp. The slip functions 
are defined by 

fc 2 ($ + *) = -8TrGa 2 pE Km D GR Si, (99) 
H- l k 2 {® + *) = -8TrGa 2 pZ Kl {f GR - 1)D G r*£100) 

in which 5, is the primordial matter density perturbation 
(measured from the CMB), -Dqr is the growth factor in 
GR, S m = D GR 5i is the GR density contrast and /gr = 
d In D G n/d In a is the GR growth rate. In the symmetron 
model (5 m = D m 6i and so 



£«m = D m / £>gr 



and 



(fm - l)*W(/< 



GR 



1) 



(101) 



(102) 
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where we have defined 



/ m = d\nD m /d\a.a = dln5 m /dhia. (103) 

We also define the linear, bias corrected, growth rate for 
the galaxies. This is defined by assuming that 



D b s tA + A 



(104) 



where Si is the initial Gaussian perturbation, is the 
bias corrected growth factor for galaxies, and the zero 
mode Ao is the source of the bias. We then define the 
bias corrected galaxy density contrast as 



cbc 



& g )S g 



where the linear bias is defined by 



(105) 



(106) 



Ao and hence the linear bias by m can be estimated directly 
from galaxy surveys using higher order statistics. We also 
define the bias corrected galaxy growth rate as 



, dln5* c 

£DC _ g_ 

' /gal ~ din a 



(107) 



In the absence of any deviations from GR, = /gr ~ 
f2° n 545 , and T, Km = S re / = 1. However in the symmetron 
model, the symmetron fifth force could result in devia- 
tions from these GR results. 

Finally, we consider the modified gravity sensitivity 
parameter defined by 



Eg = 



fc 2 (^ + $) 
-3F o 2 a-i0 gal ; 



where ga i = —5 ga \/H = —dS ga \/d\na, and so 6> ga i = 
-f^iDtA and with a = 1 today: 



Eg = 



In GR, 4 GR) - fi m0 // GR 



J gal gal 

~ O o-O- 545 

~ s 'mO"m 



(108) 



The quantities / ga i, S Km , £ K j and Eg could give us 
a feeling about how much the symmetron predictions 
deviate from the predictions made in General Relativ- 
ity. We find that these quantities depend on the red- 
shift, z, the dimcnsionless spatial scale parameter q = 
Ikl/A^f^o-ffo and the model parameters a* and z*. 
Recall that is the redshift at which the symmetron me- 
diated fifth force turns on, and a* is the relative strength 
of this extra force to that of gravity at late times (i.e., 
(1 + z) <C (1 + £*)) and on large scales. 



C. Numerical results 

Figs. 5 and 6 show / gal , £ rem , E re j and E G with z+ = 2 
and z+ = 5 respectively. In both figures we have taken 
a* = 1 and f2 m0 = 0.27. Local tests of gravity only 
constrain A2 > 10 5 -10 6 which sets the spatial scale below 



which the modification of gravity is apparent. It is clear 
from the figures that there are regions of the allowed 
parameter space where the modifications to gravity are 
potentially large. 

However, with A 2 f=s 10 6 , modification to struc- 
ture formation occur on and below Mpc scales, and on 
these scales the linear perturbation approximation breaks 
down at late times. As with the environmentally de- 
pendent dilaton model [14, 15], we expect the onset of 
nonlinear structure to suppress the fifth force giving a 
weaker observational signal, and this is confirmed in [26] . 
However, unlike the dilaton model, the redshift at which 
gravity deviates from GR, . IS cL free parameter and if 
this occurs early enough then structure will still be lin- 
ear on Mpc scales and the linearised approximation will 
hold and non-negligible changes to structure formation 
could occur. We will consider possible observational con- 
straints in the future. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have studied the background cos- 
mology and the evolution of linear perturbations for the 
symmetron model. The model utilises the so-called sym- 
metron mechanism to suppress the fifth force mediated 
by the scalar symmetron field in regions with high mat- 
ter density to evade the local gravity constraints. In the 
regions with low matter density, the fifth force is unsup- 
pressed, which could have significant consequences on the 
cosmic structure formation. 

The symmetron mechanism works because the effec- 
tive potential V B s (</>)■> which governs the evolution of <ft, 
depends explicitly on the local matter density. If the 
matter density is high, the potential has a global mini- 
mum at = about which it is symmetric; then because 
the coupling strength of the symmetron field to matter is 
proportional to <f>, the coupling vanishes and so does the 
fifth force. When the matter density drops below some 
critical value /?*, the symmetry of the potential is bro- 
ken and the minimum of V e f{ moves aways from = 0; 
because the scalar field tries to track that minimum, the 
coupling becomes nonzero and so does the fifth force. 
When inhomogeneous matter distribution is considered, 
it is possible that the symmetry in V e g is broken in a re- 
gion, but later the matter density in that region exceeds 

again because of the structure formation process, and 
the symmetry is restored so that the fifth force vanishes 
again. Obviously such a phenomenon is quite nontrivial 
to study. 

We have studied the behaviour of the symmetron field 
in background cosmology, when the symmetry breaking 
happens in the matter dominated era. We have found 
that immediately after the symmetry breaking, the sym- 
metron field lags behind its effective potential minimum 
and therefore has a negative mass-squared. This implies 
that it will experience a tachyonic period during which 
both the field and its perturbation grow exponentially. 
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FIG. 5. Effect of the symmetron field on large scale structure formation. From top to bottom the plots show: the predicted 
deviation of the (linear bias corrected) galaxy growth, /^j, from its GR value, /gr, the predicted value of the slip parameter, 
E Km , extrapolated from weak lensing measurements, the predicted slip parameter extrapolated from ISW measurements, Ei m , 
and the relative deviation of the modified gravity parameter, Eg, from its GR value. These plots are for [a* — 1, z+ = 2], and 
show the values of /g^j, E Km and Ei m at the present day, z = (solid blue line), z = 1 (dashed red line) and z = 2 (dotted 
black line) for different values of the inverse spatial scale, k. 



However this period is fairly short and the symmetron 
field will then quickly settle down to the (moving) mini- 
mum of V e g, where it oscillates. We find analytical solu- 
tions to these different phases and show that they agree 
quite well with numerical results. 

We have then investigated the general behaviour of 
the linear perturbations in the symmetron model. The 
symmetron perturbations have an effective mass-squared 
k 2 /a 2 + m 2 (4>), in which k is the wavenumber of the per- 
turbation and m{4>) is the mass of the background sym- 
metron field. On very small scales, k 2 /a 2 + m 2 (<fi) > 
even if m 2 < 0, and there is no tachyonic instability. In 
this case we notice that the growth of matter perturba- 
tion is governed by an enhanced Newton constant which 
is G = (1 + ch)Gn when the size of the perturbation is 
smaller than the Compton radius of the symmetron field 
(which means that the perturbation is within the reach 
of the fifth force), while it is the same as in GR if the 
perturbation is bigger than the Compton radius. 

On larger scales it is possible that k 2 /a 2 + m 2 (<j>) < 
given that m 2 < and the symmetron perturbation 



will experience the tachyonic growth, which naively could 
be a problem as it would mean that everything blows 
up. We find, however, that the tachyonic behaviour does 
not affect the growth of matter perturbation significantly, 
at least on sub-horizon scales. The reason is that the 
tachyonic period is so short that the even faster growth 
of <j> itself has negligible effect in cosmology, let alone 
that of S(f>. This conclusion has been confirmed by our 
numerical results. 

Finally, we have studied the structure formation in the 
linear regime under the adiabatic approximation, i.e., as- 
suming that the scalar field always tracks the minimum of 
its effective potential. Our result shows that, depending 
on the redshift z+ of the symmetry breaking, the effect 
of the fifth force on the large-scale structure could be 
important. If z+ is big, then the symmetry breaking hap- 
pens at earlier times when the matter perturbations are 
well described by the linear theory and our analysis gives 
accurate result. In such cases we find that the fifth force 
comes into effect fairly early and we expect bigger devi- 
ations from the ACDM model predictions. If z+ is small, 
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FIG. 6. Effect of the symmetron field on large scale structure formation. From top to bottom the plots show: the predicted 
deviation of the (linear bias corrected) galaxy growth, /^j, from its GR value, /gr, the predicted value of the slip parameter, 
E Km , extrapolated from weak lensing measurements, the predicted slip parameter extrapolated from ISW measurements, £i m , 
and the relative deviation of the modified gravity parameter, Eg, from its GR value. These plots are for [a* — 1, z+ = 5], and 
show the values of /i^i, S Km and Ei m at the present day, z = (solid blue line), z = 1 (dashed red line) and z = 2 (dotted 
black line) for different values of the inverse spatial scale, k. 



then the fifth force takes effect quite late, when matter 
perturbations have already entered the nonlinear regime. 
In such cases the linear theory is no longer accurate and 
better results can only be obtained using nonlinear nu- 
merical simulations [26]. 

To summarise, the symmetron mechanism is shown to 
work very successfully to suppress the fifth force to an un- 
detectable level in our solar system, while still allow non- 
negligible effects to be found in the cosmological struc- 
ture formation process. This is very important for the 
study of dark energy, because it means that one has to 
look at the small scales (relevant for galaxies and galaxy 
clusters) to find the implications of a theory which is 
supposed to modify cosmology on the very largest scales. 
Full constraint of the symmetron theory parameters will 
be studied in the future. 



Notes added: After this work had been completed, we 
noticed Ref. [24] on arXiv which also contained an inves- 
tigation of the background cosmology in the symmetron 
model. Our work mainly focus on the linear perturbation 
growth in the model, and is thus different from Ref. [24] . 
A separate work, Ref. [26], concentrates on the struc- 
ture formation of the symmetron model in the nonlinear 
regime. 
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